Preprint typeset in JHEP style - HYPER VERSION 



Jet evolution in the A/" = 4 SYM plasma at strong coupling 



Y. Hatta 

Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Ibaraki 
305-8571, Japan 



E-mail: tiatta@het . ph .tsukuba. ac . jp 



E. lancu 

Institut de Physique Theorique de Saclay, F-91191 Gif-sur-Yvette, France 



E-mail: Edmond.Iancu@cea.fr 



A. H. Mueller 

Department of Physics, Columbia University, New York, NY 10027, USA 



E-mail: 



amhOphys . Columbia. edu 



Abstract: Within the framework of the AdS/CFT correspondence, we study the time evolution 
of an energetic 7£-current propagating through a finite temperature, strongly coupled, J\f = 4 
SYM plasma and propose a physical picture for our results. In this picture, the current splits 
into a pair of massless partons, which then evolve via successive branchings, in such a way 
that energy is quasi-democratically divided among the products of a branching. We point out 
a duality between the transverse size of the partonic system produced through branching and 
the radial distance traveled by the dual Maxwell wave in the AdS geometry. For a time-like 
current, the branching occurs already in the vacuum, where it gives rise to a system of low- 
momentum partons isotropically distributed in the transverse plane. But at finite temperature, 
the branching mechanism is modified by the medium, in that the rate for parton splitting 
is enhanced by the transfer of transverse momentum from the partons to the plasma. This 
mechanism, which controls the parton energy loss, is sensitive to the energy density in the 
plasma, but not to the details of the thermal state. We compute the lifetime of the current for 
various kinematical regimes and provide physical interpretations for other, related, quantities, 
so like the meson screening length, the drag force, or the trailing string, that were previously 
computed via AdS/CFT techniques. 
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1. Introduction 

Motivated by some experimental results from the heavy ion program at RHIC, which suggest 
that the deconfined, 'quark-gluon', matter produced in the early stages of an ultrarelativistic 
nucleus-nucleus collision might be strongly interacting (see, e.g., the review papers |l], |j and 
references therein), there was recently an abondance of applications of the AdS/CFT correspon- 
dence to problems involving a strongly-coupled gauge plasma at finite temperature and/or finite 
quark density (for a recent review see [|| ) . While the early such applications have focused on the 
long-range and large-time properties of the plasma, so like hydro dynamical flow and transport 
coefficients, more recent studies have been also concerned with the response of the plasma to 
a 'hard probe' — an energetic 'quark' or 'current' which probes the plasma on very short time 
and distance scales, much shorter than the thermal scale 1/T (with T being the temperature). 
Although the relevance of such applications to actual hard probes in QCD is perhaps not so 
clear (since, by asymptotic freedom, QCD should be weakly coupled on such short space-time 
separations), the results that have been obtained in this way are conceptually interesting, in 
that they shed light on a new physical regime — that of a gauge theory with strong interactions 
— which for long time precluded all first-principles theoretical investigations other than lattice 
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gauge theory. For these results to be accompanied by conceptual clarifications, a physical inter- 
pretation for them is strongly needed, but this seems to be difficult without direct calculations 
in the original gauge theory. A different strategy, which is less rigorous, is to propose a physical 
picture based on general arguments and then demonstrate that this picture is consistent with 
all the available results (to the extent that comparisons are possible). This will be our strategy 
in this paper. 

The physical problem that we shall consider — the propagation of an 'electromagnetic' 
current (actually, an 7£-current) through the plasma — is particularly well suited for our pur- 
poses since, first, it has a strong overlap with several other problems previously considered in 
the literature and, second, it does not require any extension of the AdS/CFT correspondence 
(like the introduction of additional D7-branes) beyond the 'minimal' framework of its original 
formulation |||, [|, ||] — so it avoids any potential artifact due to such extensions. Besides, this 
problem has another useful feature: for the gauge theory of interest — namely, the conformally 
invariant N = 4 supersymmetric Yang-Mills (SYM) theory — the current-current correlator in 
the vacuum (or 'vacuum polarization tensor') is protected by supersymmetry ]7|, in such a way 
that the full result in the strong 't Hooft coupling limit A = g 2 N c — » oo is exactly the same as 
the corresponding result in lowest-order perturbation theory (i.e., the one-loop approximation). 
This property will allow us to recognize — via a comparison between the space-time picture of 
the one-loop process in the gauge theory and the 'supergravity' picture of the dual AdS/CFT 
calculation — an interesting 'duality' between the physical transverse size of the partonic system 
into which the current is evolving and the radial distance in AdS§. In turn, this duality will be 
a key ingredient of our proposal for a physical interpretation. 

Specifically, we shall follow the time-evolution of an 7£-current which in the plasma rest 
frame propagates like a plane-wave in the z direction, with a large longitudinal momentum 
q S> T and frequency u ~ q, and which at t = has zero transverse size. At weak coupling, 
this current would develop a partonic fluctuation involving just two partons — massless fields 
from M = 4 SYM which carry 7^-charge and transform in the adjoint representation of the 
color group SU(iV c ). For a space-like current, this pair would grow in transverse space up to 
a maximal size L ~ 1/Q and live for a relatively long time t c ~ q/Q 2 , thus acting like a 'color 
dipole' which can mediate the current interactions with an external target. (Here Q 2 = q 2 — u> 2 
is the virtuality of the current, and is positive in the space-like case.) For a time-like current 
(Q 2 < 0), the pair can dissociate after a time t ~ t c and thus give rise to two free partons which 
separate from each other, so like the quark and antiquark jets produced in e + e~ annihilation. 
This perturbative, one-loop, picture formally applies also to the full vacuum polarization tensor 
at strong coupling, because of the non-renormalization property alluded to above. This explains 
the polyvalence of the 7£-current as a 'hard probe', including at strong coupling 1 : by tuning the 
virtuality and the momentum of the current, we can mimic a color dipole, a 'meson', or a pair 
of jets with the desired values for the system size and rapidity, and then study the interactions 
between this partonic system and the plasma (or any other target). 

But, of course, all that applies to a current propagating through the vacuum, and there is a 
priori no guarantee that a similar strategy should also work in the thermal bath. In a previous 

1 We should also mention here the use of the current in the calculation of the rate for 'photons' or 'dileptons' 
production in the M = 4 SYM plasma at strong coupling H] . 
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analysis ||, we have considered the case of a space-like current — i.e., the problem of deep 
inelastic scattering off the plasma — and shown (via the appropriate AdS/CFT calculation) 
that this strategy still works, but only for not too high energies. A current with relatively low 
energy, such that q < Q s /T 2 , propagates through the plasma essentially without interacting, 
so like a 'small meson' in the approach of Refs. [10, 11, 12, 13, 14, 15], where the 'meson' was 
made with two 'heavy quarks' attached to a D7-brane. Note that the above condition on q can 
be rewritten as a lower bound on the current virtuality, Q 3> Q s , with Q s ~ (gT 2 ) 1//3 the plasma 
saturation momentum, and also as an upper bound on the transverse size of the effective dipole, 
L <C l/Qs, which is then consistent with the respective bound (the 'meson screening length') 



found in Refs. gg, y, |TJ. 



But for higher energies q > Q /T , or, equivalently, lower virtualities Q < Q s , the analysis 
of Ref. Q shows that the current is very rapidly absorbed into the plasma, over a time t s ~ 
q/Q 2 s {q) oc g 1 / 3 which is much shorter than the period t c ~ q/Q 2 required for the formation of 
a nearly on-shell partonic fluctuation. That is, for such a high energy, the current cannot be 
assimilated with a color dipole anymore (not even over a finite period of time), as it disappears 
before a dipole can form. The above result on t s can be restated by saying that the current 
propagates through the plasma over a longitudinal distance z s oc q 1 ^ before it disappears. 
Interestingly, the same parametric estimate has been very recently found for the penetration 
length of an effective gluon [16]. In view of the physical picture that we shall develop, this 
similarity is not just a coincidence, but rather it reflects the universality of the dissipation 
mechanism in the strongly-coupled plasma. 

In Ref. H|, the current lifetime t s has been inferred from physical considerations, but no 
temporal evolution (on top of the usual phase e~ tujt from the definition of the plane-wave) 
was explicitly considered. In this paper we shall extend that analysis by addressing the time- 
dependent problem, for both space-like and time-like currents, with the purpose of elucidating 
the dynamics of the current and the mechanism responsible for its dissipation. As already 
mentioned, we shall consider initial conditions such that the current has zero transverse size 
at t = 0. In the dual string calculation, this is represented by a vector field propagating in 
the ^dSs-Schwarzschild metric which at t = is localized at the Minkowski boundary r — > oo 
of AdS§. For t > 0, this perturbation propagates inside the bulk of the '5th dimension' as a 
Maxwell wave, i.e., according to the Maxwell equations in the yldS'5-Schwarzschild geometry. 
For the kinematics of interest, these equations can be suggestively rewritten in the form of 
time-dependent Schrodinger equations, to be presented in Sect. g. 

We first analyze the zero-temperature case (in Sect. |^), where the relevant geometry is 
purely AdS§. By solving the 'Schrodinger equations' in the approximations of interest and 
comparing the results to the space-time picture of the quantum fluctuation of the current into a 
pair of massless fields, we find that the two pictures match indeed with each other after provided 
one identifies the radial dimension r in AdS§ with the inverse 1/L of the transverse size of the 
partonic fluctuation: more precisely, R 2 /r «-> L, wit R the curvature radius for AdS$. For 
instance, on the supergravity side, the temporal scale t c ~ q/Q 2 appears as the time after which 
the wave has penetrated into the bulk up to a distance r ~ R 2 Q. This corresponds, on the gauge 
theory side, to the fact that the partonic fluctuation requires a formation time t c ~ q/Q 2 in order 
to grow up to a transverse size L ~ 1/Q. Notice that this formation time is a genuine quantum 



- 3 - 



effect, which reflects the uncertainty principle. This suggests that, for the problem at hand, 
the 5th dimension of AdS$ somehow mimics, within the context of the classical supergravity 
calculation, the phase-space for physical quantum fluctuations in transverse space. 

This r *-* L duality helps us identify a physical picture for the current dynamics in the 
vacuum. The naive one-loop picture (2-parton fluctuation) cannot be right at strong coupling, 
where any of the two partons produced in the original splitting of the current can further radiate, 
or scatter off the vacuum fluctuations. Quantum field are known for their strong propension to 
radiation. When the coupling is weak, the radiation is concentrated in specific corners of the 
phase space (collinear and soft radiation for a gauge theory), where the smallness of the coupling 
is compensated by the phase-space available to radiation. This leads, e.g, to the celebrated 
DGLAP and BFKL evolution equations in perturbative QCD. But at strong coupling, there is 
no reason why radiation should be restricted to small corners of phase-space; this should most 
naturally proceed via the 'democratical' branching of the original quanta into two daughter 
partons with more or less equal shares of their parent energy and momentum. We therefore 
propose a physical picture for current (or parton) evolution a strong coupling (and in the vacuum) 
in terms of successive branchings, such that the energy and the virtuality are divided by two (for 
simplicity) at each individual branching, and that the lifetime of a generation is determined by 
the formation time for the next generation, in agreement with the uncertainty principle. As we 
show in Sect. |5.2| , this simple picture reproduces indeed the results of the respective AdS/CFT 
calculation, at a qualitative level. 

Moving now to the finite-temperature case, we find (by solving the appropriate 'Schrodinger 
equations', in Sect. ||) that, for low and moderate energies, such that q -C Q 3 /T 2 , and for not too 
large values of time, the dynamics remains essentially as in the vacuum. For a space-like current, 
this confirms the previous results in Ref. Q: a 'small meson', with size L ~ survives 
almost unaltered in the plasma. (The 'meson' can decay via tunnel effect, but the decay rate is 
exponentially suppressed M.) But for a time-like current, the vacuum-like evolution proceeds 
only up to a time tf ~ ^/7/T, with 7 = q/Q 3> 1 the Lorentz factor of the current. At this time 
tf, the 'current' — or, more precisely, the partonic system produced via its successive branching 
- has propagated inside the plasma over a longitudinal distance Zf ~ yfy/T and has extended 
in transverse space up to a size Lf ~ l/^yT. Note that this value Lf is precisely of the order 



of the meson screening length computed in Refs. |1Q, [L2j], and more commonly written as 
Lf ~ (1 - t> 2 ) 1 / 4 /T, with v the velocity of the meson. More detailed comparisons between our 



approach and previous results in the literature will be performed in Sect. 5.1 



For larger times t >tf, this partonic system starts interacting with the plasma, although its 
transverse size is still much smaller than 1/T. In this regime, the dual string calculation shows 
that the Maxwell wave feels the attraction of the black hole and thus undergoes an accelerated 
fall towards the horizon, that it gets close to after a time of order tf. When this radial dynamics 
is translated to transverse space, via the r <-> L duality mentioned before, it implies that the 
partons feel a constant, decelerating, force in the transverse directions, 

dp± T 2 n -n 

~dT ~ ' (L1) 

whose precise physical meaning remains a little mysterious to us (see the discussion in Sect. [5^). 

By including the action of this force on the dynamics of branching, we will be able to show that 
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medium-induced branching is a likely physical scenario, which is qualitatively consistent with 
the results of the AdS/CFT calculation. By 'medium-induced branching' we mean that the 
change in virtuality over the lifetime of a partonic generation is due to the transfer of transverse 
momentum to the plasma, at a constant rate ~ T 2 (see Sect. |5.2| for more details). 

The same physical scenario — an accelerated fall of the Maxwell wave into the black hole, 
which physically corresponds to medium-induced branching — holds also in the high-energy 
regime at q > Q 3 /T 2 , for both space-like and time-like currents, but in this regime, this scenario 
starts to apply much earlier, after a time t s ~ q/Q 2 (q) which is much shorter than the formation 
time t c ~ q/Q 2 for an on-shell partonic fluctuation. (We recall that Q s ~ (qT 2 ) 1 / 3 is the plasma 
saturation momentum.) Thus, such an energetic current disappears in the plasma before having 
the time to create 'jets' (on-shell partons). Within this scenario, we find a natural physical 
explanation for the 'drag force' experienced by partons in the plasma (originally computed for a 
heavy quark probe in Refs. |l7], [l8|]), and also for the 'trailing string solution' (the string attached 
to an energetic heavy quark propagating through the plasma ^WSs-Schwarzschild geometry [17, 
|j~8| ] ) , that we recognize as the dual of the enveloping curve of the spatial distribution of partons 
produced through medium-induced branching. It remains as an interesting open problem to 
understand whether this scenario can naturally accommodate also other quantities computed 
within AdS/CFT, so like the jet quenching [19], or the transverse momentum broadening |2^, gjj . 

One can succinctly summarize the previous discussion as follows: While a highly virtual (and 
hence very small) space-like 'meson', with virtuality Q 3> Q s (q), can survive in the strongly- 
coupled plasma (essentially without feeling the latter), this is not also the case for the jets 
produced by the decay of a highly virtual time-like current, which can only propagate over 
a longitudinal distance zj ~ \fl/T before disappearing in the plasma. Less virtual (or more 
energetic) partonic systems, for which Q <C Q s (q), cannot form in the first place : the virtual 
partons melt in the plasma over a time t s ~ q/Ql(q), which is too short to get on-shell. The 
dissipation mechanism is universal, i.e., the same for all kind of partons — off-shell or on-shell, 
massive or massless, quarks, scalars, or gluons — and it consists in medium-induced branching. 

Let us conclude this Introduction with a remark concerning the spatial distribution of the 
hadrons produced via the decay of a time-like current in the vacuum. For the concept of 
hadron to make sense, we supplement the theory with an infrared cutoff A and assume that the 
current is 'hard' relative to this cutoff : \Q\ A. Then the decay of the current via successive 
branchings will produce a system of ~ IQI/A hadrons with small transverse momenta ~ A, 
which are isotropically distributed in the transverse plane. (In the rest frame of the current, this 
distribution would be spherically symmetric.) This observation is the 'time-like' counterpart of 
a previous result that, at strong coupling, there are no large-x partons (i.e., no partons carrying 
a sizeable fraction of the hadron longitudinal momentum) in the wavefunction of an energetic 
hadron |2^, ^3|, [|. Rather, all partons have fallen down at very small values of x (smaller than 
x s ~ T/Q in the case of the M = 4 SYM plasma Q and, respectively, x s ~ (A/Q) 2 for a hadron 
|23|j), via successive branchings within space-like cascades. Because of that, a high-energy 
'nucleus-nucleus' collision which would liberate these partons would produce no jets, but only a 
multitude of particles at central rapidity which are isotropically distributed in transverse space, 
with small momenta p± ~ A. This picture is similar, and possibly related, to a very recent result 
1 24 ] showing that the energy distribution of the particle produced in the strong-coupling analog 
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of the e + e~ annihilation exhibits spherical symmetry. On the other hand, this picture looks 
quite different from the corresponding one in perturbative QCD and also from the respective 
results for nucleus-nucleus collisions at RHIC. 



2. General equations 

We would like to address the following physical problem: at time t = 0, an 7£-current with 
momentum q oriented along the z axis and energy u acts in the M = 4 SYM plasma, producing 
a system of SYM quanta which then propagate through the plasma until they finally disappear. 
In the dual gravity problem, this current dynamics is represented by the propagation of a 
Maxwell-like gauge field A m in the background geometry of the AdS§ black hole (representing 
the M = 4 SYM plasma) Q . The corresponding metric reads 

ds2 = i^TBl 2 + dx2 + X_ du 2 (2>1) 

u 4ii z /(u) 

where T is the temperature of the black hole (the same as for the M = 4 SYM plasma), R is the 
curvature radius of AdS$, t and x = (x,y,z) are the time and, respectively, spatial coordinates 
of the physical Minkowski world, u is the radial coordinate on AdS$, and f(u) = 1 — u 2 . Note 
that our radial coordinate has been rescaled in such a way to be dimensionless: in terms of the 
more standard, dimensionfull, coordinate r, it reads u = (ro/r) 2 , with vq = irR 2 T. Hence, in 
our conventions, the black hole horizon lies at u = 1 and the Minkowski boundary at u = 0. 

Since the radial component A u vanishes at u = 0, it is possible and convenient to work in 
the gauge where A u is identically zero. Then, the non-zero components A^, with fi = 0, 1, 2, 3, 
are of the form: 

A^t,x,u) = e~ M « z A^u) , (2.2) 

where A^it, u) has a relatively weak time dependence, such that {dtA^l/A^ <C oj, and obeys 
the initial condition that, for t — ► 0, Au(t, u) is localized near u = 0. (The precise structure of 
the initial condition is not really needed.) The relevant equations of motion are the Maxwell 
equations in the AdS§ Schwarzschild geometry. They are most conveniently written as a set of 
equations for the transverse components Aj, with i = 1,2, and, respectively, the longitudinal 
component a(u) = d u Ao(u) (the z component A z is not independent, but it is related to Aq via 
the equations of motion), and read j|] 

4' + ^ A - Z^H. ( k 2 f +^)a i = (2.3) 

and, respectively, 



/ 1 uf 2 V J dt 2 



a"+ &la'- 4^ (k 2 f + ^)a = 0, (2.1) 



uf uf 2 \ dt 2 , 

where a prime on a field indicates a u-derivative and we have introduced the following, dimen- 
sionless, variables 

w = , k = -2— , t = 2nTt , z = 2-kTz . (2.5) 
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(The variables vj and z will be used later on.) The above equations can be further simplified for 
our present purposes. First, as already mentioned, the field Au(t,u) is slowly varying in time; 
hence we can write 



^ A, « e— ^ {-^ - ) A, , (2.6) 

where we have neglected the second-order time derivative of A^. Furthermore, as we shall see, 
the interesting dynamics happens near to the boundary at u = (and hence far away from the 
horizon at u = 1), so we can replace / — > 1 in the coefficients of the above equations everywhere 
except in the term k 2 f = k 2 (l — v?) : indeed, in that term, the small quantity u 2 <C 1 can be 
amplified by the longitudinal momentum k 2 of the incoming current, which becomes very large 
at high energy. Notice that this term k 2 u 2 = fc 2 (ro/r) 4 oc /c 2 T 4 /r 4 represents the potential for 
the long-range (in r) gravitational interaction between the current and the black hole. Since we 
keep only this particular medium effect in our equations, but neglect those which would signal 
the black-hole singularity at u = 1 (the gravity-dual hallmark of a thermal system), we expect 
our subsequent results to apply to matter distributions more general than a finite-temperature 
plasma: similar results should hold for any matter distribution which is infinite and homogeneous 
(say, a cold matter), after replacing the energy density ©oo oc N 2 T 4 of the finite-temperature 
plasma by the corresponding quantity for the cold matter. 

After these simplifications, we end up with the following equations for A\ and a, valid when 
u < 1 : 



d 

2iw-, Ai = -uA" ± K 2 Ai - k 2 u 2 Ai 

dt 

d 

2iw—zd = -ua" -a! ±K 2 a-k 2 u 2 a, (2.7) 
dt 

where we have introduced the notation K 2 = \k 2 — w 2 \ and the plus (minus) sign in front of 
K 2 corresponds to a space-like (time-like) current. Note that, with its above definition, K 2 is 
always a positive quantity. 

It will be furthermore useful (especially in view of constructing approximate solutions) to 
rewrite these equations in a form which resembles the time-dependent Schrodinger equation. 
To that aim, we shall perform the following changes of variable and functions: 

X = 2v^, a(t,u) = —^(t,x), A = ^X(f>i(t,x)- (2.8) 

Also, we shall be mostly interested in the high energy regime where u ~ q S> Q, meaning 
w ~ k 3> K, so that we can replace w ~ k in the coefficients of the equations. Then our final 
equations, valid for high energy and x<l, read 



dt V 2k d X 2 8/cx 2 2k 32 



8t V 2k dx 2 8k X 2 2k 32 
To avoid a proliferation of cases, we shall mostly focus on the longitudinal case, cf. Eq. ( |2.10| ) 
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As anticipated, we choose initial conditions such that, at t = 0, the fields are localized near 
u = (i.e., x = 0)- We shall implement that by working in Fourier space; we thus write, e.g., 



i>(t, X ) = I dee- m ^(e, X ), (2.11) 

where ip(s, x) is a wave packet in e peaked around e = 0. For convenience, we shall take this 
wave packet to be modulated by a Gaussian, that is, 

^(t, X ) = J dee- i£i -£f(e, X ), (2.12) 
where the width a is constrained by 

K 2 

k » a » — (2.13) 

k 



and the additional ^-dependence in ty(e,x) will be fixed by Eq. ( p. 10 ). The first inequality in 



Eq. ( |2. 13| ) (k 3> a) ensures that the typical energy fluctuations obey e <C w ~ k, as originally 
assumed. The second inequality (a S> K 2 /k) is necessary to allow for quasi-localized configura- 
tions in the initial condition and at early times (see Appendix A). This implies that the allowed 
fluctuations can be quite large, \w — k\ ~ K 2 /k, so that the space-like, or time-like, nature 
of the current (depending upon the sign of w — k) becomes apparent only for sufficiently large 
times, after the effects of the initial condition have dissipated. In fact, as we shall see in the 
next section, the early-time behavior of the solution ip(t,x) is independent of K 2 , and hence 
the same for both space-like and time-like currents. 

Another boundary condition refers to the behavior of the solution in the stationary regime 
(meaning, for large enough times) and at relatively large values of x '■ since a black hole is a 
purely absorptive medium, from which no signal can escape, the solution i/j(t, x) near the horizon 
located at x = 2 must be a purely outgoing wave (i.e., a wave departing from the boundary and 
impinging in the black hole). In fact, this outgoing-wave behavior should manifest itself already 
for relatively low values x ^ 1> an d hence can be used as a boundary condition on Eq. (|2.10| ), 
since the potential term in this effective Schrodinger equation is monotonous everywhere except 
near the boundary at x = 0, and hence it cannot generate reflected waves. This will be further 
discussed in the forthcoming sections. 



3. Jets in the vacuum 

We start with the zero-temperature case, i.e., with the propagation of the Abelian current 
through the vacuum. The results to be obtained here will not only serve as a level of comparison 
for the subsequent discussion of a plasma, but they will also prepare that discussion, at two levels 
(at least): First, as we shall see, there are special regimes (so like early times) where the dynamics 
in the plasma is quite similar to that in the vacuum. Second, the physical interpretation of our 
results, whose understanding is the main purpose in this paper, is easier to introduce and 
motivate in the context of the vacuum dynamics. 
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3.1 The vacuum polarization tensor 

We begin this discussion with the stationary case, i.e., the case where the fields are purely plane 
waves in the physical, 4-dimensional, space, meaning that the fields denoted with a tilde in the 
previous discussion (e.g., A^) are independent of time. The corresponding AdS/CFT calculation 
will provide the retarded current-current correlator (or 'polarization tensor') in Fourier space 
(with = (w,Q,0,q)) : 

R„u(q) = i J d 4 xe-^^(x )([J At (x),J 1 ,(0)]), (3.1) 

and in that sense it is the analog at strong coupling of computing momentum-space Fourier 
diagrams in perturbation theory at weak coupling. More precisely, R t±u (q) is obtained by differ- 
entiating the classical action with respect to boundary values of the fields at u = : 

H - (9) = a£lv (3 ' 2) 



where = A^ {u = 0) in the notations of Eq. (|2.2| ) and S is the four-dimensional action density, 
which is homogeneous: S = Jd 4 xS = AVAtS, with AV At = the volume of space-time. In 
turn, the classical action density can be fully expressed (after using the equations of motion) in 
terms of the values of the field An(u) and of its first derivative at u = : 

A " ' ~ r -A d u A* + fA 3 d u A* 3 + , (3.3) 



s 



16 



u=0 



where the z component d u A 3 is determined by the EOM as d u A% = —(w/kf)d u Ao. A star on a 
field denotes complex conjugation: the classical solutions develop an imaginary part (in spite of 
obeying equations of motion with real coefficients) because of the outgoing-wave condition at 
large u (see below). Via Eq. ( |3.2[ ), this introduces an imaginary part in R fML ,(q) which physically 
describes the dissipation of the current in the original gauge theory. In fact, the imaginary part 
of the expression within the square brackets in Eq. (3.3) is independent of u (as it can be checked 



by using the EOM) and hence it can be evaluated at any u (2? 

Note that, even in this zero-temperature context, we keep using the dimensionless variables 
u and x, which were previously defined with respect to the black-hole horizon tq = ttR 2 T. 
Throughout this section, it will be understood that ro and T are arbitrary (length and, respec- 
tively, momentum) scales of reference, which are used to construct dimensionless coordinates at 
intermediate stages of the calculations, but which will drop out from the final, physical, results. 
In this context, the radial coordinates u and x can take on all the values from to oo. 

We now turn to the actual calculation of the vacuum polarization tensor for the N = 4 SYM 
theory in the strong coupling limit. The outcome of this calculation is already known (the case 
of a space-like current has been treated in detail in Ref. ||), so our subsequent presentation 
will be very streamlined, with emphasis on the physical interpretation of the results. We shall 
give details for the longitudinal sector alone. The relevant equation reads (cf. Eq. ( |2.4| ) in which 
we let / -» 1 and d 2 /dP -» -w 2 ) 

1 K 2 

a" + - a! T a = , (3.4) 

u u 
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where we recall that the upper (lower) sign in front of K 2 = \k 2 — w 2 \ corresponds to a space-like 
(time-like) current. After the change of variable x = ^y/u, this is recognized as the equation for 
the v = Bessel functions, of either real, or imaginary, argument (depending upon the sign in 
front of K 2 ). The solution is constrained by the boundary condition (see Ref. Q for details) 

lim [ua'(u)] = k(kA + wA 3 )\ u=Q , (3.5) 

together with the condition of regularity at u — > oo. For the space-like case, these conditions 
uniquely determine the solution as 

a(u) = -2k(kA + vuA 3 )\ u=Q K (2Ky/u) (space-like) . (3.6) 

The other independent solution, Iq{2K^/u), is rejected since it would exponentially diverge as 
u — ► oo. For the time-like case, on the other hand, the general solution is a superposition of 
oscillating Bessel functions, 

a(u) = ciJ (2i^) + c 2 N (2K^) , (3.7) 

so the condition of regularity at u — > oo introduces no constraint. To fix the solution in this case, 
we shall require a(t, u) = e~ lu)t a{u) to be an outgoing wave at large u; then, the solution becomes 
imaginary, with the appropriate sign for the imaginary part to yield the retarded polarization 
tensor, via Eq. Q3.2| ). This constraint implies c\ = —ic2 which together with the boundary 
condition ( |3.5| ) completely fixes the solution as 

a(u) = -i7rk(kA + mA 3 )\ u ^H§\2K\/u) (time-like), (3.8) 

where Hq 1 ^ = Jo + «No is a Hankel function encoding the desired outgoing-wave behavior at 
large u : a(t,u) oc e - t ( ujt - K x) w hen = 2^/u » 1/K. 

The transverse-wave solutions can be similarly obtained, but the above solutions for the 
longitudinal wave are in fact sufficient to complete the calculation of R fJ/U (q): by Lorentz and 
gauge symmetry, the polarization tensor is transverse (with = (—1, 1, 1, 1)) : 

Rfj,v(q) = [r]^u ~ "qt! r (Q 2 ) (vacuum) , (3.9) 

and the scalar function R(Q 2 ) can be computed with the longitudinal waves alone. A standard 
calculation, which involves the removal of a logarithmic divergence in the real part at u = (the 
AdS/CFT analog of ultraviolet renormalization), finally yields 

R(Q 2 ) = ^||r ( ln " ^(-Q 2 H^?j , (3.10) 

where Q 2 = —uj 2 + q 2 is the current virtuality in physical units and fi is the renormalization 
scale. As expected, the imaginary part is non-zero only for a time-like (Q 2 < 0) current, which 
can decay into the massless fields (adjoint scalars and Weyl fermions) of the M = 4 SYM theory 
carrying 7£-charge. 
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As anticipated in the Introduction, this result (|3.10| ) is exactly the same as the respective 
perturbative result to one-loop order, so in particular its imaginary part (formally) describes the 
current decay into one pair of massless fields (a quark-antiquark pair, or two scalar fields). This 
interpretation is, of course, only formal: At strong coupling, the current can couple to arbitrarily 
complicated multi-particle states, but it so happens that, for the N = 4 SYM theory, the total 
cross-section is determined by the two-particle final states alone. Our present calculation being 
an 'inclusive' one — it provides the total cross-section, but it does not discriminate between the 
various final states — , its result cannot give us any direct insight into the nature of these final 
states. We shall later try to gain such an insight based on physical considerations. 

We conclude this discussion of the stationary case with another point of physical interpreta- 
tion, which points towards an interesting 'duality' between the radial dimension in AdS$ and the 
transverse size of the current (or, more precisely, of the partonic system into which the current 
has evolved) in the physical Minkowski space. Consider a space-like current, for definiteness. 
The modified Bessel function in Eq. ( |3.6[ ) decays exponentially when 2K^[u 2> 1, meaning that 
the current penetrates in the radial dimension only up to a finite distance no ~ 1/AK 2 , or 
Xo ~ 1/K (recall that x = %y/u)- The more virtual the current is, the closer it remains to the 
boundary. This is quite similar to the picture of the current in transverse space, as familiar 
in perturbation theory (for either QCD or J\f = 4 SYM): the virtual current fluctuates into 
a quark-antiquark pair (a 'color dipole') whose transverse size is inversely proportional to the 
current virtuality: L ~ 2/Q. This analogy is in fact even closer: the longitudinal wave solu- 
tion in Eq. ([O]) involves the same modified Bessel function as the wavefunction describing the 
dipole fluctuation of a longitudinal photon in lowest order perturbation theory, which allows 
us to identify the respective arguments as Kx <-> QL/2. Recalling that K = Q/2irT for a 
space-like current, we deduce the correspondence x nTL. The same argument holds in the 
transverse sector, where Ko is replaced by Kx- Also, the argument can be adapted to a time-like 
current, since the oscillatory Bessel functions in Eq. ( |3.8| ) are rapidly decaying at radial distances 
u » 1/4K 2 . 

The above writing of the correspondence between the transverse size of the partonic fluctu- 
ation and the radial coordinate in AdS§, namely, x kTL, explicitly involves the temperature 
and hence it might look a bit formal in the present context of the vacuum (but this writing will 
be natural for the subsequent discussion of a plasma). To avoid confusion, it is preferable to 
recall the definition x = 2(ro/r) with ro = itR 2 T in order to rewrite this correspondence in the 
form 2R 2 /r «-> L, from which the temperature has dropped out. 

To summarize the previous discussion, the vacuum polarization tensor in the N = 4 SYM 
theory at strong coupling formally describes the fluctuation of the current into a pair of elemen- 
tary, massless, fields, whose transverse size appears to be in a one-to-one correspondence with 
the radial distance for the current penetration in AdS$. As we shall see, this interpretation is 
comforted by the discussion of the time-dependent case, to which we now turn. 

3.2 Jet evolution in the vacuum 

As explained in Sect. |2|, we are interested in the time evolution of vector fields A^(t, u) which at 
t = start as a perturbation localized near u = (or x = 0). Given the correspondence between 
X and L (the transverse size of the current), as argued at the end of the previous subsection, we 
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see that this initial condition corresponds to a current which, at t = 0, is point -like in transverse 
space. To describe its evolution, we shall use the Schrodinger form of the equations of motion, 
cf. Eqs. ( |2.9| )-( |2~To| ), and focus on the longitudinal sector, for definiteness. The vacuum version 
of these equations is obtained by removing the last term oc x 4 111 the potential; this yields, e.g., 



dip 



i d 2 



1 



2k dx 2 8kx 2 



K 2 \ 



(3.11) 



for which we shall construct solutions obeying the relevant initial and boundary conditions. To 
remain as simple as possible, we shall consider approximate solutions which are valid piece- 
wise in x an d which are sufficient to illustrate the main points of physics. A more systematic 



method to construct solutions to Eq. (3.11), which is based on the wave-packet decomposition 



in Eqs. (2.12)-(2.13), will be described in Appendix A. 

We first note that at early times, so long as x remains smaller than a critical value Xc — 1/2K 
(the corresponding limit on time will be determined later on), the second term, proportional to 
K 2 , in the potential in Eq. ( |3.11| ) can be neglected compared to the first term cx 1/x 2 - Thus, 
this early-time dynamics is identical for both space-like and time-like currents. With the K 2 
term omitted, Eq. (3.11) admits the following, exact, solution 



Mix) 



■ Vx ih>£. 

-l^r- e 2t 
t 



This implies that the actual longitudinal wave (cf. Eq. (2.8)) 

a(t, x) 



1 1 

— — ip oc — e at 



(3.12) 



(3.13) 



is localized near x = at t = (although this is not exactly a delta function). With increasing 
time, the energy density carried by the wave fl3.13| ) diffuses towards larger values of %■> so that 
the typical distance traveled by the corresponding wave-packet after a time t is 



Xdiff(t) 



(3.14) 



This behavior holds so long as Xdis(t) < 1/2K, meaning for times t < t c ~ k/K 2 . In physical 
units, this yields t c ~ q/\Q 2 \, which is precisely the coherence time of the high-energy current; 



that is, this is the time interval which controls the Fourier transform 2 in Eq. (3.1), as it can be 
checked by rewriting the complex exponential there as 



e -iu)t+iqz ^ g - 



-iq(t-z)+iQ 2 t/2q 



(3.15) 



where we have used u — q — Q 2 /2q at high energy. 

It is interesting to compare these results with the evolution of the quark-antiquark ('dipole') 
fluctuation of a virtual photon in perturbation theory in QCD: if the photon dissociates at t = 



2 At least in the vacuum, i.e., in the absence of other time scales which are introduced by a medium. As we 
shall see in Sect, m a finite-temperature plasma introduces a new such a scale indeed. 
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into a point-like pair of fermions, then with increasing time the transverse size of this pair is 
increasing diffusively, due to quantum dynamics, like [26, 27 1 



L~y?, (3.16) 

until it reaches a maximal size L ~ l/-v/|Q 2 | at a time t c ~ q/\Q 2 \- (To avoid cumbersome 
notations, we shall often write Q instead of and Q 2 instead of \Q 2 \ within parametric 

estimates. For instance, the coherence time will be estimated as t c ~ q/Q 2 , where the modulus 
on Q 2 is implicit for a time-like current.) For t > t c , the pair is either recombining back into 
a photon, or — if the photon was time-like — it splits apart, thus giving rise to two on-shell 
particles which move away from each other. 

Clearly, the early time (t < t c ) evolution of the qq fluctuation in perturbation theory is 
very similar to the corresponding evolution of the vector perturbation in AdS§ provided one 
identifies % ~ TL, in agreement with the discussion at the end of Sect. [Oj. As we show now, 
this correspondence persists also for larger times t > t c . 

Indeed, for t 3> t c , one can heuristically estimate the time-derivative in the l.h.s. of 
Eq. (3.11) as d/dt <C l/i c ~ K 2 /k (this heuristic argument will be confirmed by the wave- 



packet analysis in Appendix A). That is, the time-derivative term in the equation is much 
smaller than the last term ~ K 2 /k in the potential, which becomes the dominant term when 
X S> 1/2K. Hence, for t 3> t c and x ^> Xc, the equation simplifies to 

d 2 

V = ±K 2 ip, (3.17) 



dx 2 

(i.e., Schrodinger equation in a flat potential) with the obvious, acceptable, solutions 

{e~ Kx (space-like) 
(3.18) 
e lKx (time-like) . 

which are time-independent and coincide, as they should, with the asymptotic versions (valid at 



large x ^ °f ^ ne respective stationary solutions constructed in Sect. 3J. In the space-like 

case, the above solution confirms that, for times t > t c , the perturbation remains localized near 
the boundary, within a distance X ^ 1 /K. In the time-like case, it implies that the actual wave 

-i-wi+ikz f -r /7 ~\ ■ ^ 



a(t, z, x) = e~ mt+lkz -- oc exp <^ -ik(t - z) - i— t + iK X > (3.19) 
Vx { 2k ) 

propagates with constant group velocity v„ = Kj k 1 along the radial direction of AdS^ : 



d (K 2 \ K ~ 



Note that, for i ~ t c , we have Xg ~ ^-/K, as it should for consistency with the previous solution 
at early times. Hence, after the wave packet has diffused up to a distance Xc ~ 1/K inside the 
bulk of AdS§, it furthermore propagates with constant radial velocity, so like a free particle. 
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This free-motion pattern at t > t c looks, of course, natural, in view of the flatness of the 
potential in Eq. ( [3.17 ). In Appendix B, we shall verify that Eq. ( |3.20| ) is the same as the radial 



part of the geodesies of a massless classical particle which propagates in AdS§ with longitudinal 
velocity v z = q/tu and radial velocity v x = Q/uj, which is the same as the group velocity in 
Eq. ( 3.2C| ) (recall that id ~ q); this classical particle is massless since v% + v^ = l. On the other 



hand, the classical particle dynamics cannot reproduce the diffusion at early stages (a genuinely 
quantum effects), nor the fall of the wave into the black hole (to be later described, in Sect. 



Via the correspondence x kTL, the result in Eq. ( 3.2C| ) is again consistent with the 



transverse dynamics expected for a time-like current which dissociates into a pair of massless 
particles (in lowest-order perturbation theory). Indeed, this result translates into 



L(t) ~ 2^ — t = 2J\-v 2 z t = 2v ± t, (3.21) 

Id 

where v z = q/uj is the common longitudinal velocity of the two massless particles, as inherited 
from the current, and v± = \J\ — v 2 is the modulus of their transverse velocity: the two particles 
move in opposite directions in the transverse plane, so the transverse distance between them 
increases like L = 2v±t. The above results have been obtained by working in the high-energy 
regime where Ld ~ q 3> Q. However, it is easy to repeat the analysis for other regimes, with 
similar conclusions. For instance, for a zero-momentum current {q = 0, Q = Ld), one finds that 
Xgif) = t, i- e -i the two particles move in the transverse plane at the speed of light (v± = 1). 

At this point, one should again stress that the reason why the physical picture looks so simple 
in the strongly-coupled N = 4 SYM theory is because of the non-renormalization property 
of the current-current correlator, as alluded to before. The AdS/CFT calculation correctly 
provides the total cross-section for the current decay, but it does not capture the detailed 
nature of the final states. Formally, the total cross-section is saturated by the two-particle 
final state; therefore, it is the dynamics of this particularly simple state which emerges, via the 
correspondence x ^ TfTL, from the dual calculation on the gravity side. 

This discussion has an interesting corollary: it shows that the current can be also viewed 
as a device for introducing a pair of elementary, massless, fields of the M = 4 SYM theory 
(Weyl fermions or adjoint scalars) at a given radial distance within AdS$, which is controlled 
by the current virtuality: x ~ 1/K or r ~ QR 2 . This is tantamount to fixing the transverse 
size L ~ 2/Q of the partonic pair in the physical space. That is, the dynamics of the current at 
times larger than the coherence time is the same as that of 'meson', or of a 'color dipole'. For 
a space-like current, this effective 'meson' simply sits at x ~ meaning that its transverse 

size is fixed. For a time-like current, this 'meson' propagates with constant velocity along the 
radial direction %, meaning that its transverse size grows at constant speed. In the next section 
we shall study the influence of a thermal bath on this dynamics. 

4. Jets in the plasma 

We now turn to the problem of main interest for us here, which is the propagation of the current 
and of the associated partonic system through a strongly-coupled N = 4 SYM plasma with 



2 
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temperature T. We are interested in 'hard probes', so we shall choose a current with relatively 
high virtuality (either space-like, or time-like) : Q = \/\Q 2 \ T (or K 3> 1), which probes 
the structure of the plasma on distances much shorter than the thermal wavelength 1/T. We 
shall mostly consider a relativistic current, for which q ~ u 3> Q, but the non-relativistic case 
(q <C uj) will be briefly discussed too, for completeness. In fact, the physically most interesting 
case — the one where the medium effects should be truly relevant — is when the coherence 
time t c ~ q/Q 2 of the current is much larger than 1/T, so that the current explores a relatively 
large longitudinal slice of the plasma, with width Az ~ t c S> 1/T. This implies a lower limit on 
the current momentum: q 3> Q 2 /T, which is tantamount to the condition that the Bjorken-x 
variable 3 , defined as x = Q 2 /2qT (in the plasma rest frame), be very small: x -C 1. 



4.1 Physical regimes 

We would like to determine the characteristic time scale for the dissipation of the current in the 
plasma; in the dual, gravity, problem, this is the time scale for the fall-off of the Maxwell field 
An into the black hole. As we shall see, this scale is controlled by the dynamics at relatively 
small u <C 1 , and thus is insensitive to the detailed geometry of the black hole near its horizon 
at u = 1. As explained in Sect. B the strength of the gravitational interactions between the 
wave and the black hole is proportional to the wave longitudinal momentum q, and also to the 
temperature. In view of that, we shall be led to distinguish between two important physical 
regimes: 

(i) a relatively low-energy (or low-temperature) regime at qT 2 <C Q 3 (or k -C K 3 ), where 
the medium effects are strongly delayed, so that the current dynamics proceeds as in the vacuum 
up to time scales much larger than t c , and 

(ii) a very high-energy (or high-temperature) regime at qT 2 3> Q 3 (or k 3> K 3 ), in which 
the current dissipates very fast, on a time scale much shorter than t c . 

One can understand these various regimes by studying the potential in the Schrodinger-like 
equations ( |2 . 9|) ( pTTo| ) . Let us first recall from the previous section that the time-dependence in 
these equations is important only at very early stages, when the typical values of x are so small 
that the only relevant term in the potential is the first term, oc l/k\ 2 , which is independent of 
both the current virtuality and the properties of the medium. Hence, for such early times, the 
equations describe diffusion, for both space-like and time-like currents, and in the same way 
as in the vacuum. But for later times, where the meaning of 'later' is generally medium- and 
energy-dependent (see below), the solutions approach a stationary regime, where Eqs. (|2.9f) — 



(2.10) can be replaced by their time-independent versions, of the generic form — ip" + Vip = Eip 
with E = (since E corresponds to the time-derivative, which is negligible). Hence, this late- 
time behavior is determined by the time-independent Schrodinger solution with zero energy in 
the potential V(x)- 

In what follows we shall study this solution in the longitudinal sector alone (the corre- 
sponding discussion for the transverse sector being very similar). The respective potential is 
graphically illustrated in Figs. [I], ||, and ||, for the various physical regimes: low energy and 



3 This variable is especially relevant for a space-like current which undergoes deep inelastic scattering off the 
plasma H. 
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Figure 1: The potential in Eq. ( 2.1(\ ) for the space-like case at relatively low energy (k <C K 3 ). 

viz) 
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free streaming 



medium branching 



Figure 2: The potential in Eq. (2.1C ) for the time-like case at relatively low energy (k <C K 3 ). 



space-like in Fig. [j], low energy and time-like in Fig. ||, and high energy (both space-like and 
time-like) in Fig. ^. Some general features of the dynamics are already clear by inspection 
of these figures: At sufficiently small values of x, the potential is the same as in the vacuum; 
the medium starts to be felt only at relatively large values of x, where the potential is a de- 
creasing function which describes attraction by the black hole. Also, in the high-energy regime, 
the dynamics is essentially the same for both space-like and time-like currents. Finally, in the 
space-like case at least, there is an obvious difference between the low-energy regime, where 
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the potential barrier constrains the wave to remain on the left of the 'classical turning point' 
at x ~ 1/K (so like in the vacuum), and the high-energy regime, where the barrier has dis- 
appeared and the wave can easily propagate towards the horizon. (The transition between the 
two regimes occurring at k ~ K 3 corresponds precisely to the situation where the height of the 
potential barrier becomes negligible.) In the time-like case, on the other hand, there is never 
such a barrier, so the difference between the low-energy and high-energy regimes looks perhaps 
less obvious. As we now explain, this difference is nevertheless important in that case too. 

To that aim, it is convenient to separate the three pieces in the potential in Eq. ( 2.10| ), 
which play different physical roles. We thus write (for a time-like current) 

1 K 2 A-v 4 

v( *' = -8i?- »--£- = * + v fl + ifc. («) 

Let xab denote the value of x at which Va and Vb become comparable with each other, and 
similarly for xac an d Xbc- We have, parametrically, 

1 1 [K , a , 

Xab ~ , xac ~ , Xbc ~ y -q • (4.2) 

Consider now the two physical regimes alluded to before: 

(i) Low energy regime : k < K 3 xab < Xac < Xbc 

It is easy to check that, in this regime, the potential admits the piecewise approximation 

Va + V B for % < Xbc 
V( X ) oc { (4.3) 
Vc for x > Xbc ■ 
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Thus, so long as the relevant values of x (those where most of the wave energy is located) remain 
much smaller than xbc-, the dynamics is the same as in the vacuum. This is the situation at 
sufficiently small times, that we shall study in the next subsection. But when x ^ XBC th e 
medium effects become important and entail the fall of the wave into the black hole. This fall 
be studied in Sect. |4.3| . Note that the transition point x ~ Xbc is where the potential changes 
from a plateau to a rapidly decreasing function (see Fig. Q). 

(ii) High energy regime : k » K 3 => xbc < XAC < XAB 
In this regime, the virtuality-dependent piece Vb is never important, and the potential can be 
approximated as (see also Fig. ||) 

{V A for x < XAC 
(4.4) 
Vc for x > XAC , 

which now holds for both space-like and time-like currents. Unlike the low-energy potential, 
the high-energy one has no intermediate plateau, but only a pronounced peak at x ~ XAC- in 
fact, the transition between the two energy regimes, which occurs for k ~ K 3 , corresponds to 
the situation where the two endpoints of the plateau, xab and xbc-, merge with each other, and 
also with the position xac of the emerging peak. 

4.2 Early time dynamics 

In this subsection, we shall consider the dynamics of the current at early times, which is insen- 
sitive to medium effects (except for its validity limits, as introduced by the plasma), and thus 
can be inferred from the corresponding discussion in Sect. ||. 

We start with a simple case, which has not been explicitly covered by the previous discussion, 
but which can be viewed as a special limit of the low-energy case: a non-relativistic, time-like, 
current with q <C u and u S> T. The corresponding 'Schrodinger' equations are obtained from 
Eqs. ( |2 . 9| ) ( pTlQ ) by omitting the last, medium-dependent, term in the potential (since this 



term is comparatively small for any < x ^ 2) and replacing k — > w and K 2 — > vo 2 in all the 
other terms. The ensuing equations are the same as in the vacuum, and therefore so is also 
the current dynamics, until it eventually dissipates into the plasma. Specifically, at very early 
times t < l/w, the associated Maxwell wave diffusively penetrates into AdS§, up to a distance 
X ~ (meaning that the current diffusively spreads out in the physical space, up to a size 
L ~ Then, for larger times t 3> 1/w, the wave propagates inside AdS§ with constant 

group velocity v g = 1, meaning that the current has decayed into a pair of massless particles 
which move away from each other at the speed of light: Xg(t) = ^ or L(t) = It. This behavior 
continues until the wave hits the horizon (xg = 2 or t ~ 1) and is thus absorbed by the black 
hole. Physically, this means that, for times t ~ the separation between the products of 
the current decay has become as large as the thermal wavelength, L ~ 1/T, so these products 
cannot be anymore distinguished from the plasma fluctuations — they get 'lost' in the plasma. 

We now turn to a relativistic current with w ~ k ^> K 2 , for which we have to distinguish 
between two physical regimes, as previously explained: 

(i) At moderate energies, such that K 2 C fc < K 3 , and for a time-like current, the upper 
limit xbc ~ y K/k for the vacuum-like behavior is much larger than xab ~ which is 
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the maximal penetration length through diffusion. Hence, in this regime, one can observe both 
stages of the evolution identified in Sect. |3|: a 'diffusive' stage at early times t < i c ~ k/K 2 , 
during which the wave-packet progresses according to Eq. ( 3.14| ), and a 'free motion' stage at 



t > t c , during which the wave has constant group velocity v g = K/k. This second stage lasts 
until the position Xg(i) = v gi of the wave-packet becomes comparable to Xbc ) this happens at 
a time tf ~ \JkjK (the subscript / on tf stands for "free motion"), which represents the upper 
time limit for this vacuum-like dynamics. In physical units, this yields a time scale 




1 

Q 2 " f 



t f ~ ~J± » tc ~ ^ » -, (4.5) 



which, as indicated above, is much larger than the coherence time. 

The physical interpretation of this result, as deduced via the correspondence x *~~* kTL, is 
as follows: the time-like current develops a partonic fluctuation with size L ~ 1/Q over a time 
of the order of the coherence time t c , and then decays into a system of massless particles which 
move freely (without feeling the plasma) up to a time t ~ tf. At this stage, the transverse extent 
of the partonic system has increased up to a value 




L f = 2v ± t f « -, (4.6) 



which is still small as compared to the thermal wavelength in the medium. For even later times 
t > tf, the dynamics is driven by the medium and will be analyzed in the next subsection. 

It is interesting to notice at this point that, although it influences the dynamics at late 
times t > tf, the plasma has no effect on the polarization tensor, which is still given by the same 
expression as in the vacuum, namely Eqs. (3.9)-( j3T0|) with Q 2 < 0. Indeed, if one considers a 



(time-like) plane-wave solution, so like in Sect. [Oj, then for x *C Xbc this solution is again given 
by Eq. fl3.8| ), which describes an outgoing wave for any x within the range xab <C x ^ Xbc- 
The change in the solution at larger x ^ Xbc has no incidence on the calculation of R fJil/ (q), 
which is determined by the behavior of the solution near x = 0, cf. Eq. (|3.2|), Physically, this is 
so since the current disappears by decaying into a pair of massless fields, so like in the vacuum, 
and the subsequent fate of these fields is irrelevant for the calculation of the total decay rate 
(the imaginary part of R^u). This is similar to the calculation of the total cross-section for 
e + e _ annihilation in lowest-order perturbation theory: this cross-section is fully given by the 
e + e~ annihilation into a qq pair, although the quark and the antiquark are not the actual final 
states in the experiments. The 'late stages' radiations or interactions involving the quark and 
the antiquark, although essential for hadronisation and the composition of the final state, do 
not affect the total e + e~ cross-section. 

(ii) At ultrarelativistic energies k 3> K 3 , the medium effects (as enhanced by the energy) 
become important already at the very short radial distance xac ~ 1/A; 1 / 3 <C 1/K, meaning that 
the current has no time to fully develop its partonic fluctuation before getting absorbed. At 
very early times, such a fluctuation starts to develop via diffusion, but this process is interrupted 
after the relatively short period t s ~ A; 1 / 3 , or 

i. ~ i (|) « U ~ (4,, 
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(which is still bigger than 1/T though), when the wave-packet has diffused up to xac- Physically, 
this means that the growth of the fluctuation is stopped at a transverse size 

L -~\/f ~ Hi) " 3 << r r (4 - 8) 

which is much smaller than the natural size ~ 1/Q for the same fluctuation in the vacuum. 
Accordingly, the dynamics is very different from the vacuum case, in the sense that the dissi- 
pation of the current and the respective polarization tensor are now controlled by the medium. 
This polarization tensor, which is now identical for space-like and time-like currents, has been 
computed in Ref . @ , and the result was used to deduce a partonic interpretation for the struc- 
ture of the plasma in its infinite momentum frame. In the next subsection, we shall address 
this problem from a different perspective, by following the time evolution of the Maxwell wave 
during its fall towards the black hole. The insight that we shall gain in this way will allow us to 
propose, in Sect. ||, a physical picture for the current dissipation in the strongly coupled plasma. 

4.3 The fall of the wave into the black hole 

We now come to the most interesting physical situation, which refers to the (relatively) late 
stages of the current evolution in the plasma and the mechanism for current dissipation. In this 
situation, the dynamics is controlled by the medium-dependent piece in the potential (the piece 



V c ~ x i n Eq. (4.1)), which yields the following equation of motion in Schrodinger form 



.fy ( 1 d 2 k x 4 * 



dt V 2k dx 2 32 
As previously explained, for a time-like current this equation holds both at moderate energies, 
K 2 <C k <C K s , and in the high-energy limit k 3> K 3 , but in ranges for x which are different 
in the two cases (cf. Eqs. ( f4.3| )-(fO|)). For a space-like current, it holds only for k 3> and 
X » XAC- 

We have not been able to find an exact solution to this equation, but we shall construct a 
WKB approximation to it. To that aim, we use the wave-packet representation in Eq. ( 2.12j ). 
The corresponding WKB solution reads then 

^(*>X) = / dee~ t£{t ~ to) ~^ — ]— exp 
J Pe{X) 




p £ ( X ) = ^/lk{e - V c (x)) , (4.10) 

where we have kept only the outgoing wave and xo is either xbc, or xac, depending upon 
the physical regime under consideration (and similarly for to). At this point, we notice that 
the interesting values of t are large enough for the typical energy e ~ 1/t to be much smaller 
than the potential \Vc\ ~ For instance, when k 3> K 3 , we are interested in t > i s ~ fc 1 / 3 
together with x ^ Xac ~ 1/A; 1//3 . Then we can expand the square root within p £ to linear order 
in e and perform the ensuing Gaussian integration, to obtain 

x 



^{t,x) ^ —r-r exp 
PoKX) 



1 ./",,, , , <y 




T 2 



(t-k) 



1 1 

Xo X 



(4.11) 
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We have also used here 



dp £ 
de 



x 



4 - u& 



2\V C \ x 2 J de 



\[---) (1.12) 

\xo x; 



Xo 



The modulus \ip(t, x)\ tells us where the energy of the wave-packet is located at time t. Clearly, 
with increasing time, the peak of the energy distribution moves along a trajectory 

*® - T^W^> ■ <" 3 > 

which can be recognized as the trajectory of a classical particle with mass m = k and zero total 
energy moving in the potential Vc(x)- This is so almost by construction (because the classical 
trajectory defines the group velocity for the WKB solution), and can be also checked by starting 
with the respective particle equation of motion, that is, 

k^ = -^ = k^, (4.14) 
dt 2 d X 8 ' V ' 



whose zero-energy solution brings us back to Eq. ( 4.13] ), as it should 



The approximation ( [4.13 ) holds only so long as \ <C 1, but it can be used to estimate the 



duration of the fall of the wave-packet in the potential; namely, x(t) reaches the horizon at 
X = 2 when t ~ th, with 

4 4 

t h -to = 2 ~ — . (4.15) 

Xo Xo 

Note that, parametrically, xo ~ l/^o f° r both the intermediate-energy and the high-energy 
regimes. Thus, % — to ~ io> so the total time th necessary for the wave-packet to propagate 
within AdS^ from the boundary to the horizon is of order to. Essentially, half of this time is 
used to (rather slowly) reach the ' point of no return' at xo *^ 1 (the point where the wave 
starts to feel the attraction of the black hole) and the other half to travel along the considerably 
larger distance from xo U P to the horizon. Clearly, for x > Xo the motion of the wave-packet 
is accelerated by the gravitational force. In fact, Eq. ( |4,14| ) implies that the group velocity 
X = dx/di approaches the speed of light (from below) when the wave gets close to the horizon: 

X 2 

X = — X — 1 when x — 2- (4-16) 

In reality, the particle cannot cross the horizon (in the laboratory frame), but only asymptotically 
approach to it, so its velocity can never become exactly one. But when x — 2, the terms in 
the potential which have been neglected in writing Eq. ( |4.S| ) become important and modify 



Eq. ( 4.16 ). Let us also note here the 'dual' version of the above equation, which is particularly 
suggestive : using the correspondence x ~ TL and introducing the parton transverse momentum 
p± via the uncertainty principle, i.e., p± ~ 1/L, we arrive at 

^~-T 2 . (4.17) 
dt 1 ; 
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That is, there is a transverse force acting on the colored partons in the plasma, and this force is 
independent of the parton momentum and of order T 2 . This result will play an important role 
in the physical discussion in Sect. 5.2. 

Let us conclude this section with a brief summary of the previous results, in physical terms: 

(i) At moderately high energies, such that Q/T <C q/Q <C (Q/T) 2 , a time-like current 
creates SYM quanta in the plasma at t = 0, which then pass through three stages before 
disappearing : 

(a) a relatively short period of diffusion, at t < t c with t c ~ q/Q 2 , during which the current 
gets replaced by a bunch of on-shell massless partons (only two such partons being manifest in 
our calculation) , with overall transverse size L ~ 1/Q ; 

(b) a longer period of free streaming, from t ~ t c up to tf ~ (l/T)^Jq/Q, during which the 
partonic system propagates without feeling the plasma and expands in transverse space up to a 
size L ~ (l/T)y/Q/q <C 1/T, and 

(c) an equally long period, during which the interactions with the plasma cause an acceler- 
ated expansion in transverse space, up to a maximal size L ~ 1/T; when this size is reached, the 
partons move apart from each other at the speed of light (v± ~ 1) and disappear in the plasma. 

The lifetime of the current is controlled by the two last stages, and thus is of order tf, which 
is much larger than the coherence time t c . 

(ii) At very high energies, such that q/Q ^> (Q/T) 2 , a virtual current, either time-like or 
space-like, goes through only two stages before dissipation: 

(a) a relatively short period of diffusion, at t < t s with t s ~ q/Q 2 (q) <C t c , during which 
the partonic fluctuation coming from the current grows up to a transverse size L ~ 1/Q s (q), 
which is however not sufficient for the partons to become on-shell, and 

(b) an equally short period of accelerated expansion, at times t > t s , which lasts until the 
partonic fluctuation reaches a transverse size L ~ 1 /T and disappears in the plasma. 

Both stages contribute on equal footing to the lifetime of the current, so this lifetime is of 
order t s . 

Above, we have introduced the saturation momentum Q s (q) ~ (qT 2 ) 1 / 3 , which is the natural 
transverse momentum scale to discuss the high-energy scattering off the plasma ||. This is 
defined as the 'critical' virtuality for the transition between the low-energy and the high-energy 
regimes: we have indeed q/T ~ (Q s (q)/T) 3 . Note that, at high energy, Q s (q) replaces Q as the 
virtuality scale which determines all the physically relevant scales for transverse size and time. 
(For instance, both t c and tf reduce to t s after the replacement Q — > Q s (q)-) Moreover, Q s (q) 
is also the critical scale for the disappearance of the potential barrier in the space-like case. 



5. Physical discussion 

This section will contain no new technical developments, but in spite of that it could be viewed 
as the main section of this paper; indeed, here is where, on the basis of the previous results, we 
shall develop our physical picture. But before we do so, it is useful to make contact with some 
previous approaches in the literature. This should make clear that the physical picture that we 
shall later develop will also shed light on some of the results of these previous approaches. 
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5.1 Relation with previous approaches 

In this subsection, we shall show that our previous results are consistent with, and shed new 
light on, previous studies in the literature |l0|, [0], [L2|, 13, 14, 15, 16, 17, [l8|], concerning the 
propagation of 'heavy quarks', or of quark-antiquark 'mesons', in the strongly coupled N = 4 
SYM plasma. 

To be specific, let us start with an example (we shall present some others later on): Consider 
a time-like current at moderately high energy (Q/T <C q/Q <C (Q/T) 2 ) and focus on our 
estimate ([0|) for Lf — the transverse size at which the partonic system created via the decay of 
the current starts to lose energy in the plasma. After introducing v z = q/oj and 7 = l/yl — v% = 
q/Q, this result can be rewritten as 

which turns out to be the same as the parametric estimate for the 'meson screening length' 
^max(v z ,T) computed in Refs. [1C, 11, 12]. The latter is the maximal transverse size that a 
'heavy meson' (a bound system of a quark and an antiquark) can have in the plasma, when it 
propagates along the z axis with a constant velocity v z < 1. So long as the transverse size t of 
the meson is smaller than f max (v z ,T), the qq pair is bound indeed and it moves freely through 
the plasma — that is, it moves at constant speed v z without the need for an external force. 
Mathematically, this is described by a string connecting the quark and the antiquark through 
AdS§, which rigidly moves together with its endpoints. For £ > £ max (v z , T), on the other hand, 
this connecting-string solution ceases to exist and is replaced by two disjoint strings trailing 
behind the quark and, respectively, the antiquark, which extends in AdS§ all the way up the 
horizon. This trailing strings exerce a drag force on their fermionic sources, meaning that energy 
is transferred from the quark and the antiquark to the black hole |L7], |l8|, |2{| . 

This meson dynamics is similar to that emerging from the present calculation, except for 
the fact that, here, the transverse size of the partonic pair is not a free parameter any longer, 
but is rather fixed by the dynamics (and hence it depends upon time). Namely, over the time 
interval t c < t < tf, which is large as compared to the formation time t c ~ q/Q 2 of the pair, the 
time-like current can be effectively viewed as a 'meson' made with a pair of on-shell, massless, 
partons (Weyl fermions or adjoint scalars), which propagate together along the z axis with 
constant velocity v z and at the same time separate from each other in transverse directions 
with a (small) relative velocity 2v±, with v± = a/1 — v 2 <C 1. Hence, the transverse size of the 
'meson' increases like L = 2v±t. When this size reaches the critical value in Eq. ( [4.6| ) or ( |5.1| ) 
(this happens at a time t ~ tf ~ ^/7/T), the partons start to lose energy to the plasma and 
then the 'meson' disappears, so like in the corresponding calculation in Refs. 11, 12|. 

If this analogy is correct, then we should also find some similarity between the trailing-string 
solutions in Refs. [17, 18, 28, 29| and the dynamics of the current at late times t 3> tf, where 
dissipation is important. We will later show that this is indeed the case, but before that let us 
emphasize some differences between our actual problem and those in the previous literature: 

• Our 'partons' are massless fields from the Lagrangian of the N = 4 SYM theory. A pair 
of such fields is produced by the current, at a distance \ ~ ^/K inside AdS$ which is controlled 
by the current virtuality. Physically, this means that this pair is produced with a transverse 
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size L ~ 1/Q- By contrast, in previous studies, the 'partons' are 'heavy quark probes', that 
is, massive fermions which do not belong to J\f = 4 SYM and whose introduction requires an 
extension of the original AdS/CFT correspondence: the dual partners of these fermions are open 
strings ending on a D7 brane which extends along the radial direction of AdS§ up to a distance x 
inversely proportional to the fermion mass. This makes it difficult to consider massless, or even 
light, quarks, since the corresponding D7 brane would enter the black hole. Moreover, difficulties 
appear even for heavy quarks in the high-energy limit (v z —►!.), since the latter appears not to 
commute with the limit in which the ultraviolet cutoff is sent to infinity [19, [30], 15, ^jj ]. This 
particular problem does not appear in our formalism, where physics is completely smooth in the 
limit v x — ► 1, even in crossing the lightcone. 

• The partonic pair produced by the decay of a time-like current is, strictly speaking, not a 
'meson' (i.e., a bound state which would be stable in the absence of the plasma), but rather a pair 
of jets, similar to the qq jets created via e + e~ annihilation in perturbative QCD. However, the 
relative velocity of these jets in transverse space is much smaller than their common longitudinal 
velocity (v± = yT — u| <C v z ), so the fields remain close to each other for a long time ~ ^y/T; 
this may explain why their dynamics is so similar to that of a meson. As we shall shortly argue, 
the partonic fluctuation of a space-like (relativistic) current is even closer to an actual meson. 

• In our approach, there is no explicit string dual counterpart of the 'partons', or 'jets' (so 
like the endpoints of the string on the D7 brane in the 'meson' problem). The partons appear 
only in the physical interpretation on the gauge theory side, where they are supposed to be 
generated via current fluctuations. Related to that, the 't Hooft coupling A, which in the other 
approaches is introduced by the Nambu-Goto string action, is not present in our formalism, 
which in fact corresponds to the strict strong-coupling limit A — > oo. This feature sometimes 
complicates the comparison between our results and previous approaches. 

In spite of such differences, striking similarities persist (like the one already discussed in 
relation with Eq. ( |5.1D ), which support the idea that the two types of problems (the 'current' 
and the 'meson') are closely related. In what follows, we shall give some more examples in that 
sense, which suggest that our wave-packet propagating through AdS§ — the gravity dual of the 
plasma perturbation by a current — is in fact rather similar to the string solutions describing a 
meson, or heavy quark, perturbation in the previous approaches. 

Consider first a time-like current at moderate energies and not so late times (t < tf), 
corresponding to a relatively small meson. The dynamics of the connecting-string solution on 



the 'meson' side is characterized by two important scales in AdS$ (see, e.g., |12|, |15|]) : the 
position XtipOO of the peak of the string (the maximal penetration of the string in AdS$) and 
the radial coordinate Xv = beyond which the string action would become imaginary. For 

a relatively small meson I <C £max.{v z ,T), one finds Xtip(^) ~ T£, which is precisely the relation 
between the distance traveled in AdS§ and the transverse size of the partonic system that we 
advocated at the end of Sect. 3.2. This suggests that one can identify Xtip(^) with the position 
of our wave-packet in AdS$. Furthermore, recalling that 7 = q/Q, one sees that Xv ~ \/Q/l is 
the same as our 'point of no return' XBCi cf. Eqs. fl4,2|) and (|4.3|) : beyond this point, the wave 
cannot escape the attraction by the black hole. 

Similar identifications can be made also for a space-like current with moderate energy (q -C 
Q 3 /T 2 ): Then, 

Xtip(^) and Xv correspond to the first and, respectively, second 'classical turning 
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points' for the potential barrier in Fig. [j], i-e., XAB ~ T/Q and, respectively, xbc ~ \jQll- The 
partonic system has now a fixed size L ~ 1/Q (since the Maxwell wave gets stuck at x ^ Xab), 
and in that sense it looks even closer to a real meson. 

Returning to the case of a moderate-energy time-like current, let us now consider the 
dynamics at later times, t > tf, where we would like to make contact with the trailing-string 
solution in Refs. [17, |18|]. Notice first that the criterion for the onset of dissipation, which 



determines the 'critical' length (|5.1|), is the same in our problem and in the 'meson' problem: 
this is the condition that the distance traveled in AdS§ get close to the 'point of no return', 
which amounts to Xg ~ Xbc (cf. the discussion above Eq. (|4~5|) ) in our context, and to Xtip ~ Xv 
in the approach of Refs. ||i~0| , O , jl^j . There is an interesting difference though: whereas in the 
usual meson problem, this change of regime is abrupt and characterized by a sharp critical value 
^ma,x{v z ,T), in our approach it is much more gradual , as it is driven by the competition between 
the various terms in the potential. For x ^ Xbc (and hence t 3> tf), our solution describes 
a wave-packet falling into the black hole, along the trajectory in Eq. ( 4.13| ). This should be 



compared to the trajectory of the energy flow towards the horizon in the corresponding trailing- 



string solution. The latter can be inferred from the results in Refs. [17, [18| : the rate dE/dt 
for energy flow towards the horizon and the energy density d.E/dr per unit length along the 
dragging string have been there computed as (see 5 , e.g., Eqs. (3.19a) and (3.20a) in Ref. 

dE = ttV\T 2 v 2 z dE = V\ 

dt 2 v / l _ ^f ' dr 2 7 r v / r^fi? 2 ' 

(Recall that x = 2(ro/r) with tq = ttR 2 T, and that we consider a situation where % <C 2, i.e., 
r S> tq.) These expressions involve the 't Hooft coupling A, which is inherent in the 'heavy 
quark' problem (it enters via the string tension, which provides the overall normalization for the 
Nambu-Goto action, and also for the string energy-momentum tensor) and seems to prevent 
any comparison to our present results (which correspond to the limit A — > oo). However, the 
coupling constant disappears in the ratio 

which determines the trajectory of the energy flow in AdS^, and is the right result to be compared 



to our respective trajectory in Sect. |4.3| . And, indeed, Eq. (|5.3|), where for the present purposes 
we can take v z ~ 1, is identical with our respective expression ( [4.16 ), as it can be recognized 
after a change of variables. 

In fact, the identification between our wave solution and the dragging string can be made 
even more precise: the string world-sheet, as expressed by the function z(t, r) = v z t — £(r) (the 
string, which has a shape z = £(r) in the variables z and r, moves at constant speed v z , solidary 
with the heavy quark which is pulled by an external force) is exactly the same as the surface of 
stationary phase for our Maxwell wave solution a(t, z,r), as computed in Ref. p. Specifically, 



4 Similarly, for a space-like current, the transition between the moderate-energy and the high-energy regimes, 
which happens when xab ~ Xbc (i.e., for q ~ Q z /T 2 ), is smoothed out by the tunnel effect |5J. 
5 The authors of Ref. H 

use the notation u for the radial coordinate in AdS$, which is however not the same 



as our variable u in Eq. (2.1); rather, their definition for this variable is u = r/R 



2 
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in Ref. || we have studied the time-independent version of the above Eq. (4J3), but with a 
more general form for the potential, valid for all the values of x on the right of the 'point of 
no return', up to x = 2 (or r = ro). The plane-wave solution has been obtained in the WKB 
approximation as (see Eq. (D.4) in Appendix D of Ref. B) 



a(t,z,r) 



-iwt+iqz 



f q / 1 r -\- tq ( T \ \ ] 

exp < i — — - In h arctan — > . (5.4) 

P \27rrV2 r-r \r J J J 1 ; 



(For r ^> ro and after a change of variables, this expression reduces indeed to the complex 
exponential part of Eq. ( 4.1 1| ) . ) The condition of constant phase determines a hypersurface in 



AdS§ which can be written as z — (co/q)t = Q(r) + const, where uj/q ~ 1 is the phase velocity 
and C(r) is the same function as the shape of the trailing string, cf. Eq. (3.17) of Ref. [jl] 



or Eq. (10) of Ref. [fig]. Note the logarithmic singularity in Eq. (5.4) at r — > ro : this reflects 



the fact that, as mentioned in Sect. 4.2, the wave-packet can never cross the horizon, but only 
asymptotically approach to it, according to a law r — ro oc exp(— ^irTt). In the next subsection, 
we shall propose a physical interpretation for this profile £(r) of the trailing string. 

At this point, we should also emphasize an important difference between our present set-up 
and the 'meson' problem in the previous literature: In the latter, the transverse size t of the 
meson and its longitudinal velocity v z are independent parameters, so it is in principle possible to 
choose a velocity v z arbitrarily close to one (although in practice this meets with the problem of 



the UV cutoff, as alluded to before [19, 30, 15, |3l[]). In our approach, the 'meson' is dynamically 
generated as a fluctuation of the virtual current, which requires a formation time t ~ t c . This 
introduces an upper limit on the longitudinal momentum q up to which a nearly on-shell 'meson' 
can form in the plasma: when the 'meson' lifetime tt ~ (1/T)^J q/Q becomes comparable to 
its formation time t c — this happens for q ~ Q s /T 2 — , the partonic fluctuation melts in the 
plasma before having the time to become on-shell. Hence, a 'meson' can form only so far as 
q <C Q^/T 2 , which in turn introduces an upper limit on the longitudinal velocity v z of the meson 
thus generated: 7 <C (Q/T) 2 , or 1 — v 2 > (T/Q) 4 . But, of course, there is no corresponding 
limit on the energy of the current: when q 3> Q 3 /T 2 , the 'point of no return' is reached 
after the very short time t s ~ q/Q 2 (q) *C t c , and then the current disappears in the plasma 
before forming a meson. The late-time dynamics which is responsible for this disappearance, 
as studied in Sect. [O, is exactly the same for on-shell or off-shell partonic fluctuations : it is 



always characterized by the group velocity (|4.16[) for the falling wave-packet and by the surface 



of stationary phase surface determined from Eq. (5.4). This points out towards the universality 
of the dissipation mechanism in a plasma at strong coupling, which is the same for a colorful 
(heavy or massless) quark as for a (sufficiently energetic) colorless meson or current, although in 
the AdS/CFT calculation this dynamics may embrace different mathematical descriptions (e.g., 
a trailing string, or a Maxwell-like wave-packet). 

In the next, final, subsection we shall try to elucidate the actual physical mechanism which 
is responsible for this dissipation in the strongly-coupled gauge theory. 

5.2 Parton branching at strong coupling 

In our physical discussion so far, we have privileged the interpretation in which the time-like 
current first fluctuates (over a time of the order of the coherence time t c ~ q/Q 2 ) into a pair 
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of massless partons which subsequently undergo free motion, along straightline trajectories. In 
the vacuum, this scenario would hold for ever, and independently of the energy of the current, 
which only determines the velocities of the particles in the pair. In the plasma, this would 
hold only for moderately high energies, q <C Q 5 /T 2 , and for values of t smaller than the time 
tf ~ (l/T)y/q/Q at which the pair is large enough (in transverse direction) to feel the plasma. 
But even for larger times t > tf, the partons would still follow classical trajectories, which now 
describe (in the dual theory) the fall of the Maxwell wave into the black hole. At higher energies, 
q S> Q 3 /T 2 , this accelerated fall would show up right away after the early period of diffusion. 

Although consistent with our previous calculations and also quite appealing due to its sim- 
plicity, this physical picture cannot be fully right, for several reasons: First, at strong coupling, 
there is no reason why the current should couple to two-parton final states alone, or why these 
partons should undergo free motion, not even in the vacuum. Rather, the two partons produced 
in the first, 'electromagnetic', splitting can in turn radiate other partons via strong, 'color', 
interactions, and also strongly couple to the fluctuations of the vacuum, thus giving rise to a 
complicated, multi-partonic, final state. Second, the fact that the partons appear to follow 
classical trajectories in the physical (longitudinal and transverse) space contradicts the quan- 
tum nature of the M = 4 SYM theory. Although, in practice, we solve classical equations of 
motion in the dual gravity theory, the results thus obtained must somehow reflect the quantum 
nature of the original problem in gauge theory, that should fully reveal itself at strong coupling. 
This is already manifest in the dynamics at early stages, where our complex exponential in 
Eqs. ( 3.12 )— ( 3.13| ) can be recognized to describe, via the correspondence x kTL, quantum 
diffusion (i.e., quantum Brownian motion) in the physical transverse space. There should be a 
corresponding quantum dynamics hidden in the classical wave solution at later times t > t c , and 
in what follows we shall try to unveil this dynamics. 

This seems a priori difficult because of the non-renormalization property of the polarization 
tensor in the M = 4 SYM theory, as discussed in Sect. 3.1. It is precisely this property which on 
one hand makes the corresponding result to look so simple at strong coupling (cf. Eq. (|3.10p ), 
but on the other hand forbids any direct access to the detailed nature of the final state — and 
hence to the actual fate of the partons produced by the current. Here, we shall simply try to 
reconstitute this fate from the dual wave solution via physical considerations, notably, by using 
the uncertainty principle. 

Note first that by using the uncertainty principle alone one can correctly estimate the 
formation time for the partonic fluctuation, i.e., the coherence time t c alluded to above. Indeed, 
in the rest frame of the current, where its 4-momentum reads q'^ = (u/, 0,0,0) with u' = Q, 
the uncertainty principle requires a fluctuation time t' c > 1/Q. This becomes t c = jt' c ~ q/Q 2 
after a Lorentz boost to the frame in which the current is relativistic. 

We now turn to later times t > t c (but lower than tf in the case of the plasma), where the 
behavior in Eqs. fl3.19| )-( |3,20 ) applies, which in turn implies that the partonic system expands 
in transverse space according to Eq. ( |3.21 ). Previously, we have interpreted these results as 



describing the free motion of two classical particles, but now we shall argue that they are also 
consistent with quantum branching at strong coupling, a dynamics which is much more plausible 
given the circonstances. In this interpretation, the expansion of the packet in transverse space 
is the result of the degradation of the partons transverse momentum via successive branching, 
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in conformity with the uncertainty principle. 

Specifically, let us describe the branching process in terms of generations — a kind of mean 
field picture which should be reasonable at large N c — and start with a single parton (with 
longitudinal momentum q and virtuality Q 2 ) at to = 0. We shall consider only a 1 — » 2 splitting 
vertex, which is representative for the dynamics in a gauge theory. Also, we shall assume that 
each of the two daughter particles produced in a splitting takes a way, roughly, half of the 
momentum and virtuality of the parent particle. This assumption is natural at strong coupling, 
since there there is no need to look for special corners in phase-space, so like collinear or 
soft emissions, to enhance the probability for splitting. Thus, in the nth generation, as obtained 
after n > 1 successive branchings, there will be 2 n particles, each of them carrying a longitudinal 
momentum q n ~ q/2 n and a virtuality 6 Q n ~ Q/2 n . Since the virtuality was relatively small to 
start with (q ^> Q) , and it further decreases via the splitting, the particles in each generation are 
nearly on-shell. At strong coupling, a splitting occurs as fast as permitted by the uncertainty 
principle, so the lifetime of the (n — l)-th generation — the time t n — t n -i necessary to go from 
the n — 1 to n branching — can be estimated from the energy imbalance at the splitting vertex. 
This yields t n — i n _i ~ q n /Q n ~ 2 n (g/Q 2 ), which rapidly grows with n. The time evolution of 
the virtuality can therefore be estimated as 

Qn ~ Qn-1 Q n2 dQ(t) Q 2 (t) 

— Q n => —7— , (5-5) 

tn ~ t n -i q dt 7 

where 7 = q/Q is the Lorentz factor for the original parton. Introducing (by virtue of the 
uncertainty principle, once again) the transverse size L n ~ 1/Q n of the assemble of partons 
making up the nth generation, and similarly L(t) ~ 1/Q(t), we finally obtain 

^r = 7' (5 - 6) 

with C a number of order one whose precise value is not under control. This is consistent, 
as anticipated, with the law ( |3.21 ) obtained from the AdS^ calculation via the correspondence 



X ttTL. Incidentally, Eq. ( |5.5[ ) implies Q(t) ~ j/t, valid for t > t c - 

In the vacuum of the conformal N = 4 SYM theory, this successive branching would hold for 
ever, down to smaller and smaller values of the longitudinal momentum fraction z n = q n /qn-i 
of the produced partons. If, in order to mimic confinement, the theory is supplemented with 
an infrared cutoff A (say, in the form of a cutoff r m - m = AR 2 — implying a maximal value 
Xmax °c 1/ A — on the radial distance in the dual string theory) , then the splitting will continue 
until the virtuality of the softest produced partons will become of the order of this cutoff: 
Qn ~ A . The total duration of the decay process will then by controlled by the lifetime of the 
last generation, tjy — t^-i ~ 2 N (q/Q 2 ) ~ 7/A, where we have used 2^ = Q/A and q/Q = 7. 
The final partons produced in this process are relatively numerous (2^ = Q/A > 1) and have 
small transverse momenta q± ~ Qn ~ A, so they will be isotropically distributed in transverse 
space, within a disk with area ~ 1/A 2 around the longitudinal axis. 

But in the case of a plasma, this vacuum-like branching continues only up to a time tf ~ 
y/y/T, when the partonic system has expanded in transverse space up to a size Lj ~ 1/y/jT ^ 



"As before, the modulus on Q and on Q n is implicitly understood for time-like partons. 
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1/T. It is easy to understand this 'critical' value in the current scenario: at t ~ tf, the 
softest partons in the cascade have a typical virtuality Qn ~ VAf ~ ^ \/l an d a longitudinal 
momentum q N = q/2 N ~ 7 3 / 2 T. (H ere iV denotes the number of generations up to a time 
tf, hence, 2 N ~ QLf.) These values satisfy ~ QnT 2 , which is precisely the condition for 
the partons in this A^th generation to start interacting with the plasma. In Sect. [|, we have 
associated this condition with a highly energetic current, which starts to feel the plasma already 
at early stages during its evolution, before having the time to decay into on-shell partons. Here, 
we have so far considered a low energy current (qT 2 <C Q 3 ), which at the beginning decays in 
the same way as in the vacuum, but whose evolution yields a system of partons for which the 
condition Q^ N ~ QnT 2 is eventually satisfied (since the successive branchings lead to a faster 
decrease in Q n as compared to q ri ). The physical meaning of this condition can perhaps be 
better appreciated by noticing that the rate for the change in virtuality in this iVth generation 
(i.e., at time t ~ tf) is of the order 



dQ{t) 
dt 



T 2 , (5.7) 



t=t f 



which looks like a natural order of magnitude for the transverse force exerted by a plasma on 
colored partons. At a first sight, this might look consistent with a 'quasi-particle' picture for 
the strongly-coupled plasma, in which the quasi-particles (thermal excitations) have typical 
momenta ~ T and are randomly distributed in space, with a typical interparticle separation 
~ 1/T. But then the random scattering between the parton and these quasi-particles would 
increase the dispersion < p\> —p 2 ^ in its transverse momentum, rather than uniformly decrease 



its average momentum p± ~ Q. So, most likely, Eq. (5/7) does not describe thermal scattering, 
and its precise physical origin remains to be clarified. One should also notice that the partons in 
this iVth generation have a relatively long lifetime ~ y^y/T, so they can lose a significant part 
of their transverse momentum \Qn+i — Qn\ ~ T-J^j (before further decaying) even though the 
respective dissipation rate is rather small, |dQ/dt| ~ T 2 . 

For larger times t > tf, the dynamics proceeds in the same way as it would have proceeded 
at all times t > if the current was sufficiently energetic (q > Q 3 /T 2 ) to start with: namely, the 
partons keep branching, but this branching is now accelerated by their momentum loss towards 
the plasma, at an average rate ~ T 2 . Indeed, as already noticed in relation with Eq. (4.17), the 



AdS / CFT result ( ^4 . 1 6| ) for the fall of the Maxwell wave into the back hole is tantamount (after 
the identification \ i — y TL ~ T/Q(t)) to a decelerating transverse force 

~ -T\ (5.8) 

dt 

acting on the partons. Because of this force, the partons can now dissociate much faster than 
in the vacuum : the lifetime At n ~ q n /Qn of a parton from the nth generation (with n > N) 
is determined by the condition that, during this time, the parton be able to lose a transverse 
momentum of order Q n , at a constant rate ~ T 2 . This condition determines the parton virtuality 
as Q n ~ ((/nT 2 ) 1 / 3 , which is recognized as the 'saturation momentum' introduced in Sect. |I| 
- here evaluated for a parton with momentum q n ~ q/2 n . Note that this value for Q n is 
considerably larger than Q/2 n (the corresponding value in the vacuum after the same number 



-29- 



of generations), showing that the partons which are involved in this medium-induced branching 
are highly off-shell. 

It is also interesting to evaluate the rate dq/dt for longitudinal momentum loss in the plasma, 
i.e., the longitudinal force acting on the parton. We can write 

" Ql — ^#=-(,rr 3 , (5.9) 



t n -i n _i Qn/Ql "" dt 



which is recognized as the 'dynamical' version of the drag force computed in Refs. |17], 18 1. 
Namely, in those papers, one has considered a heavy quark moving through the plasma at a 
constant speed v z and found that, in order to compensate for its energy loss, one needs to drag 
this quark with a constant force dp/dt = (l/v z )(dE/dt) that can be read off the first equation 



( |5.2| ). Therefore 7 dp/dt oc 7T 2 , where 7 = 1/^/1 — v 2 is time-independent in the context of 
Refs. [|T7| , |i"8[| . Our force ( |5,9| ) can be rewritten in a formally similar way by introducing the 
Lorentz factor for the (virtual) parton as j(t) = q(t)/Q(t) — which however is now time- 
dependent, because of the uncompensated energy loss — and recalling that, for this parton, 
Q(t) = Q s (q(t)) ~ (qT 2 ) 1 /' 3 , so that we indeed have (qT 2 ) 2 / 3 ~ j(t)T 2 . 

Eq. ( |5.9[ ) tells us that a parton (or current) will totally lose its energy over a time interval 
^loss ~ (t/T)(qN /T) 1 / 3 , with q^ the value of its longitudinal momentum at the time where that 
parton (or current) has started to feel the plasma. For an 7£-current with moderately high energy 
(qT 2 <C Q 3 ), we have seen that qm ~ J 3 / 2 T, which in turn implies t\ oss ~ ^y/T, in agreement 
with the respective AdS/CFT result in Eq. (|4.5| ). But for a very energetic current (qT 2 3> Q 3 ), 
and also for a colored parton with any energy, q^ is the same as its original momentum at t = 0, 
and then we find ti oss ~ (1/T)(q/T) 1 / 3 , in agreement with the corresponding result in Eq. 



and also with a very recent result in Ref. [16]. 



We shall conclude this discussion with a physical interpretation for the trailing string solu- 
tion originally obtained in Refs. |p"7| , |l8| (and which, we recall, also emerges from our Maxwell 
wave solution, cf. Eq. ( |5.4| )). Namely, we shall argue that the function £(r) which defines the 
shape of the string is 'dual', via the identification r/R 2 <-> 1/L, to the enveloping curve of the 
parton distribution created via branching in the medium. This curve can be best visualized 



if one considers a stationary situation similar to that in Refs. [17, [lq| , namely the situation 
in which some external force is continuously giving energy to the system, in such a way that 
the leading particle propagates at constant speed v z < 1. The 'leading particle' is either our 
7£-current in the high-energy regime at qT 2 3> Q 3 , or some colored parton like the heavy quark 
considered in jlT], 18]. This particle creates a parton cascade via successive branchings and the 



average (gross) properties of this cascade are independent of time. The parton distribution being 
roughly isotropic in transverse space, as argued before, we can focus on a particular plane, say, 
(x,z) (see Fig. Q). The enveloping curve of the distribution is then the function Az(L) which 
relates the longitudinal separation Az between a parton in this plane and the leading particle to 
the transverse width L of the cascade at the position of that parton. To construct this function, 
it is convenient to return to the description of the parton cascade in terms of generations; this 



7 We take v z ~ 1, as appropriate for comparing to our present calculations. 
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Figure 4: The parton cascade generated by the current via medium-induced branching in the physical 
Minkowski space (represented here as the boundary of AdS$ at \ — ) and the trailing string attached to 
the leading particle (represented for \ «1 ). The latter is 'dual' to the enveloping curve of the former. 



implies 

n 

Az(L n ) ~ ^(l-v^Atj, (5.10) 

j'=i 

where L n ~ l/Q n and where Qj, Atj, and Vj denote, respectively, the typical virtuality, lifetime, 
and velocity of a parton from the jth generation, with j < n. According to the previous 
discussion, we have Qj ~ (^T 2 ) 1 / 3 , Atj ~ qj/Qj, and 

Vj ~ -S- = , ~ 1-$,. (5.11) 

Therefore, 

A,(L n) ~ f 7 - r ^ - r t t2l3 ~ t2l « ' (5 - i2) 

where we have used dj ~ dLj/Lj (since d(lnLj) ~ — (l/3)d(ln qj) ~ dj ; we recall that qi = 
g/2 J ). We thus have found that Az(L) ~ T 2 L 3 . After replacing L — > i? 2 /r within this result 
(in order to compare with the shape £(r) of the trailing string), we finally deduce Az(r) ~ 
T 2 (i? 6 /r 3 ). As anticipated, this is parametrically the same as the expansion of the function £(r) 
for r > ro (where our present calculation is supposed to apply). 
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By inverting the function above, we find that L grows with Az as L ~ (Az/T 2 ) 1 / 3 . This is 
correct so long as L remains smaller than 1/T, and hence for values of Az which are themselves 
smaller than 1/T. But the transverse width of the partonic distribution is roughly limited 8 to 
L < 1/T, since partons with momenta Q <T disappear in the thermal bath. This argument 
suggests that at larger distances Az > 1/T the parton distribution should approach a cylinder 
with diameter L ~ 1/T (see Fig. This is again similar to the corresponding behavior of the 
trailing string, which for ( > 1/T approaches asymptotically the horizon at r = tq |0, 18]. 



This 'duality' between the shape of the trailing string and the enveloping curve of the parton 
distribution implies a similar duality between our Maxwell wave in AdS§ and the wave which 
would describe the physical parton distribution in Minkowski space. Namely, the latter can be 



obtained by replacing r — > 2R /L in Eq. (|5.4|) , thus yielding 



a(t,z,L) e -*w^Nff*-H«A*(£) 



**(L) EE ^ 



1 Lp + L (L q 
— in h arctan — 

2 L - L \L 



L« = ^T- (5-13) 



4 

The phase of this wave is stationary along a paraboloid in Minkowski space which represents the 
enveloping surface of the parton distribution. The fact that this distribution can be characterized 
by its enveloping curve alone (i.e., that is does not depend separately upon the transverse coor- 
dinates x and y, with x,y < L, but only upon the points on a circle with radius L) can be simply 
understood in physical terms: all the points in a transverse cross-section of this paraboloid at 
a fixed value Az = z — (oj/q)t correspond to partons which belong to a same generation. Such 
partons move in phase with each other, as they have the same velocity. Accordingly, the phase 
of the wave must be the same for all such points, as indeed happens for the wave in Eq. ( 5.13j ). 
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A. Wave— packet evolution in the vacuum 

In this Appendix we shall consider the initial value problem for the Schrodinger-like equation 
( 3.11 ) in a slightly more systematic way, by using the Ansatz ( [2.12 ) in order to search for 



8 Actually, as pointed out in Refs. |3^, |33| the energy and momentum transfer from the leading particle 
can also induce longer range (L 3> 1/T) perturbations into the plasma, so like sound waves or Mach cones. But in 
order to see such perturbations, one has to take into account the feedback of the trailing string, or of the Maxwell 
wave, on the ^IdSs-Schwarzschild geometry, something that goes beyond the purpose of this work. 
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appropriate solutions. We recall that we are interested in solutions ip{t,x) which for t are 
localized near x = 0. The precise form of the initial condition is however irrelevant, and in what 
follows it will be more convenient (since mathematically simpler) not to stick to a particular 
initial condition, but rather identify a particular solution and then verify that this solution 
satisfies indeed the condition of localization at early times (t — > 0) . 

We consider the time-like case, i.e., we take the minus sign in front of K 2 in Eq. ( 3.11| ), 



since this is the most interesting case. After inserting the Ansatz ( 2.12| ) into Eq. ( p. 11 ), one can 



recognize the Bessel equation for *$>(e, x)> U P to some trivial changes of function and variable. 
An interesting particular solution is 

i>(t, X ) = j dee- t£i -^^x~,Mxv / K 2 + 2ke), (A.l) 

where the reason for choosing the particular Bessel function Jo (rather than the most general 
combination C1J0 + C2N0) should become clear in a moment. We shall estimate the above integral 
separately for small and, respectively, large values for the external coordinates t and x — the 
change in behavior occurring at t ~ i c ~ 2k / K 2 (the coherence time) and for x ~ Xc — 1/2K 
(the position of the wave-packet at t = t c ). In this analysis, one should keep in mind that the 
width a of the packet obeys to Eq. ( |2,13|) . 

Specifically, for small times t <C t c (with t > 1/a though) and for x ^ Xc the integral is 
dominated by energies e such that K 2 /2k <C e <C a, as we shall shortly check. Then one can 
neglect the Gaussian factor inside the integrand and also the K 2 term inside the argument of 
the Bessel function. The resulting integral can be performed exactly (see formula 6.631-6 in Ref. 
|35|])) and reproduces the corresponding solution found in Sect. (cf. Eq. ( |3.12| )) : 

d £e -^VxJo(xV2fc) = —pe^T . (A.2) 

In fact, it was precisely this possibility, to exactly perform the above integration, which motivated 
our choice for the Bessel function Jq in Eq. (|A.1|). 



As discussed in Sect. 3.2, this solution implies that the energy diffuses along a trajectory 
X ~ ytjk. For values of x along, or near, this trajectory, we have xV2ke ~ Vte, and then the 
integral in Eq. (A.2) is controlled by e ~ 1/t (since this is the only scale inside the integrand). 



Hence, for i <^ i c , we have e ^> K 2 /2k, as anticipated. Since, moreover, the time variable is 
restricted to t > 1/<t, the condition that e <C a is satisfied as well. 

In the other interesting case, namely for relatively large times and radial coordinates, t > t c 



and x ^ Xc the integral in Eq. (AJ) is dominated by relatively small values of e, such that 
2ke <C K 2 . Moreover, in this regime we have Kx S> l, so one can use the asymptotic expression 
for the Bessel function, to obtain 

where we have kept only the outgoing-wave component of the solution. (The other component 
would yield a solution propagating from the bulk towards the boundary.) The dominant contri- 
bution can now be obtained by expanding the square root in the exponent to linear order in e 
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and then performing a Gaussian integral : 

' dse- iei -^ +iKx+i ^ E ~ e iK x-4(*-*xf . (A.4) 



The strength of the wave, which shows where the energy is located, has a peak at t = {k/K)x- 
Hence, this solution represents a wave-packet with propagates in AdS^ with constant group 
velocity : 

K 1 



X = v g t, v g = — = - = a/1 - v 2 z , (A.5) 

k 7 



in agreement with the simpler analysis in Sect. 3.2. 
B. Classical particle falling in AdS 

In this appendix we consider the motion of a classical pointlike particle falling down into the 
AdS 5 black hole and compare its motion to the propagation of a time-like current, as previously 
studied in Sect. ||[ (See Ref. |36|] for a similar study.) 
A massless particle propagates along a null geodesies 

= ds 2 oc - ( 1 - ^] dt 2 + dx 2 + dy 2 + dz 2 + ^ , (B.l) 

V XhJ ( 27r T) 2 (l-4) 

where xh = 2. For a particle with velocity (the variables q and to are introduced below to 
facilitate the comparison with the discussion in Sect. 



dz _ q 
dt LO 



(B.2) 



in the z-direction, ( p.l| ) reads 

In the static case, v = 0, one can integrate explicitly 

t = tanh" 1 — + tan" 1 — . (B.4) 

XH XH 

The particle reaches near the horizon x ~ Xh hi time t ~ 1/T. This is indeed similar to the 
evolution of a time-like current with q = 0. 
When g ^ 0, ( jB.3[) can be rewritten as 



dt) A xlJU 2 x 4 h 



When 1 > Q/u> » X /Xhi 01 



X«2 A /^~ X , (B.6) 
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one easily finds 



dx Q 

r^j 

dt ~ LO ' 

which is equivalent to Eq. ( |3.20| ). This is valid until a time 

t < — J — ~ tf , 



(B.7) 



(B. 



when the 'point of no return' x ~ Xu is reached. The above estimates for and Xv are indeed 
consistent with the corresponding ones derived from the current dynamics in Sect. ||| (where Xv 
was rather denoted as xbc)- Therefore, the intermediate stage of the evolution of a time-like, 
moderate energy (q <C Q 3 /T 2 ), 'photon' looks indeed like a massless pointlike particle falling 
into the bulk. On the other hand, the initial, diffusive regime of the photon evolution has no 
classical analog. Moreover, a pointlike particle stalls in the x direction at x ~ X/> as clear 
by inspection of Eq. ( |B.5[ ), whereas a photon (string) keeps on falling further into the bulk as 
discussed in the main text. 
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